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Abstract 

Temperature dependence of the upper critical field H C 2 (T) for the superconducting 
magnesium diboride, MgB 2 , is studied in the vicinity of T c by using a two-band Ginzburg- 
Landau (G-L) theory. The temperature dependence of Hc 2 (T) near T c exhibits a positive 
curvature. In addition, the calculated temperature dependence and its higher order derivatives 
are also shown to be in a good agreement with the experimental data. In analogy with the 
multi-band character of Eliashberg microscopic theory, the positive curvature of H c2 (T) is 
described reasonably by solving the two -band of G-L theory. 
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I. INTRODUCTION 

Recent discovery of superconductivity in an intermetallic compound of magnesium 
diboride, MgB 2 , has initiated an intensive research activity (see for example Refs.[l-5]) in the 
condensed matter community, simply because the compound was of interest to the researches 
in many years. In addition to its record high critical temperature at around 40 K for a binary 
compound, its isotope effect was also of interest. The isotope effect was immediately studied 
by Bud'ko et al.[6] to test the microscopic theory of superconductivity. As a superconductor, 
the mechanism of superconductivity in MgB 2 involves interesting phenomena like two- 
phonon scattering and multiple gaps (Ref.[5] and the references therein), giant anharmonicity 
and nonlinear electron-phonon (e-ph) coupling [2] . 

MgB 2 is of a type II superconductor with the G-L parameter K=26 as given in Ref. [7]. 
Meanwhile, low values of critical fields H c i and H c2 are reported and the estimated values for 
H c i and H c2 at absolute zero are about a few hundreds Oe [8] and twenty thousand Oe [9-10], 
respectively. Large positive curvature of the upper critical magnetic field H c2 near critical 
temperature T c was reported. The similar positive curvature has also been observed 
experimentally in nonmagnetic borocarbides of Y(Lu)Ni 2 B 2 C [11,12]. Microscopic 
mechanisms of superconductivity in MgB 2 and discussions of the superconducting properties 
implicitly involve the employment of Isotropic Single Band (ISB) model [13]. However, 
when analysing the superconductivity using the ISB model, the multi-band character and 
anisotropic Fermi surface of the superconductor have been ignored. This leads to an 
inadequate analysis. Although the superconductivity in MgB 2 is phonon mediated, it is likely 
that an analysis beyond the single band Eliashberg theory is needed. Hence, two-band 
Eliashberg theory [14] can be complementary to the ISB model. Recently two-band 
Eliashberg model was applied [15] for determining the behaviour of H c2 in MgB 2 . 
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Although extensive theoretical studies at microscopic level have been carried out 
immediately after the discovery of MgB 2 , it is necessary to provide additional information on 
its superconducting properties by using macroscopic theory of G-L. In this paper, we propose 
that the experimental data of H c2 for MgB 2 and its differentiation with respect to temperature 
can be described in the vicinity of T c by using a two-band macroscopic theory of G-L in a 
similar way used for nonmagnetic borocarbides in Ref.[16]. 



II- THEORY 

We write the G-L free energy functional for two-band superconductors with two 
coupled superconducting order parameters of the following form similarly as in [17], 

F sc = Jd 3 r(F 1 +F 12 +F 2 + H 2 /87t) (1) 
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Here m i denotes the the effective mass of the carriers belonging to the band i (i=l,2). F; is 
the free energy of the separate bands. Fi 2 is the interaction energy term between the bands. 
The coefficient a depends linearly on the temperature T as 0Cj =Yi(T-T c j), while the 

coefficient (3 is independent of temperature, y is the proportionality constant. The quantities 8 
and 8i describe the interband mixing of two order parameters and their gradients, respectively. 
H is the external magnetic field , $ is the magnetic flux quantum and A is the vector 

i i2 

potential. In Eqns. (2) and (3), it is assumed that the order parameter y¥A are slowly varying 



functions in real space. By minimization of the free energy of Eqn.(l) 



SF/9¥*=0 , 8F/&¥*=0 (4) 
we then obtain the usual basic equation for the description of the two-band superconductivity. 
For simplicity, we assume a vector potential A = (0,Hx,0) in one dimension, then 

linearisation of Eqn.(4) gives 

2 2 2 2 2 
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4m 2 dx 2 l 2 dx 2 l 2 

Where l 2 = (hc/2ell) is the so called magnetic length. The coefficients 8i and 8 are 

determined by the microscopic nature of the interaction between different band electrons. 
Their signs can be taken arbitrarily. If the interband interaction is ignored, the Eqns. (5) and 
(6) are decoupled into two ordinary G-L equations with two different critical temperatures. In 
general, regardsless of the sign of 8, the superconducting phase transition occurs at a well 
defined temperature exceeding both T c i and T c2 (see Eqn.(8)). Without linearisation, the G- 
L equations would include a cubic term. For analytical solutions, one can make an 



approximation by replacing the cubic term of Eqn.(4) by an expression of *Pj = 



vj/2 
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To our best of knowledge, two band G-L equations have not yet been solved by taking into 
account the average of the order parameter. 

The upper critical field H c2 can be obtained as the lowest eigenvalue of the equation 
describing two-coupled harmonic oscillators in a similar way as for a single band G-L theory 
[18]. Hence, we look for eigenfunctions for the equations (5,6) of the form: 

¥j =X il exp(-ax 2 /2)+X i2 exp(-bx 2 /2) (7) 



Substituting these solutions into the Eqns. (5,6), we finally obtain an equation involving the 
upper critical field H c2 implicitly as: 
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The ratio of order parameters at T c can be found from the following equations: 
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Using Eq.(8) one obtains for the normalized upper critical field h c2 - H c2 / H c2 (°) 
h c2 = (e - c + (AG 2 - B0 + c 2 , ) J 1 2 )ao 1 • (10) 



Where = (1-T/T C ), H c2 (0) = cT c (Yjmj + y 2 m 2 )/#e . The curvature of the upper critical 



field with the dimensionless parameters reads 
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The parameters defined in Eqns. (12-15) are all dimensionless quantities. Near T c , we have an 
asymptotical behaviour for h c2 of the form: 
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Ill- RESULTS AND DISCUSSIONS 

The experimental data for the upper critical field h c2 of MgB 2 can be described with 
high accuracy by the simple expressions 
tt O l+oc 

h c2 =4^ = — 2 r witho)=(i-0)0 1+ot (i7) 

H c2 (0) (l-(l + a)co+^Gr +mco 3 )' 
This formula of Eqn.(17) was apphed for fitting the experimental data of h c2 in nonmagnetic 
borocarbides by Drechsler et al. [19]. They found an unusual positive curvature near T c . 
Miiller et al [10] have only ignored the denominator of Eqn.(17) and only used the remaining 

simple formula, 1+oc for fitting their experimental data for MgB 2 . The critical exponent a 
determines the temperature dependence more effectively. However, in such an 
approximation, they did not attempt to explain the negative curvature at low temperatures. 
Meanwhile, the saturation (negative curvature) at low temperatures is well described by the 
ratio of lira , which is assumed in analogy as sensitive to the electronic structure. For £ > m 
the temperature region with negative curvature becomes wide. For £ < m we have an upper 
critical field changing almost linearly with temperature. 

In figure 1, we plot the experimental data of Bud'ko et al. [9] with open circle symbol 
and Miiller et al. [10] with filled circle symbol for H c2 (T)in MgB 2 versus the reduced 

temperature T/T c and also show the theoretical fits to this data by using Eqns.(lO) and (17). 

The best fit to the experimental data was obtained by using Eqn.(17) (shown by solid lines) 

* 

with fitting parameters |1qH c2 (0) = 16.2 T , a = 0.23, £ = 3,m = -1. In the figure, the dotted 

line display the two band G-L fitting by using Eqn.(10). For this fitting, G-L parameters 
together with other fitting parameters are as I 

A = 0.59,B = -0.11,c =0.20,x = 3,|l H c2 (0) =15.55 T,T cl =20 K,T c2 =10K. 
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At low temperatures, H c2 is changing almost linearly with temperature. As T c is approached 
from below, linearity begins to disappear giving a shallow like dependence. As shown in 
Fig.l, there exists a reasonable agreement between theory and experiment. Note that the 
single band G-L theory would only give a linear temperature dependence for H c2 . The small 
deviations between the data and the theory may suggest some additonal effective bands to be 
considered in the G-L theory to look for the exact temperature dependence of H c2 . For 
Miiller's data, we again observe the similar temperature dependence. 

In figure 2, we plot the first and second order derivatives of h c2 of Eqn.(17) and 
Eqn.(ll) with respect to temperature versus reduced temperature T/Tc. Two band G-L theory 
gives a negative and an almost unchanging the first derivative below T c . However, as T c is 
aproached from below, the first derivative increases until a finite value at T c . On the other 
hand, the first derivative of the fitting formula Eqn. (17) starts from zero at absolute zero, 
goes through a negative region and returns to zero at T c . There is a good agreement in the 
temperature range from 20% to 95% of T c . The plots of the second derivatives are also given 
in the same figure. A good agreement in the same temperature range is observed for the 
second derivatives, too. Note that the G-L theory holds only at temperatures close to T c . 
Hence, we concentrate on the upper side of the curves rather than the lower side. 

In the case of no intergradients of order parameter (r\ =0), the G-L curvature reaches 
its maximum at the point = B/2A ~ 0.5. It should be noted that a similar result (in the 
framework of two band BCS theory) was obtained in [20], but only for a rather high value of 
the mass ratio parameters x=20. The inclusion of a negative intergradient interaction shifts 
this maximum to the region = 0. Physically, it means that in the vicinity of T c , when 
both order parameters are small , their interaction becomes crucial for the behaviour of the 
upper critical field. As shown by Eqn. (9) , the critical temperature of around T c =40 K can be 
obtained from T c i=20 K and T c2 =10 K approximately with the interaction parameter 8 = 3/8. 
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Upper critical field h c2 is governed by the parameter e l . Our fitting corresponds to T) = - 
0.16. We also take into account that the average Fermi velocity over hole tubular band is 
v F1 =4.5 10 7 cm/s, while for the remaining last Fermi surface sheets we can choose v F2 =11. 
10 7 cm/s . In the case of nonmagnetic borocarbides, we have the following parameters : v F1 
=0.85 10 7 cm/s, v F2 =3.8 10 7 cm/s [14]. The theoretical fit to the experimental data of 
nonmagnetic borocarbides by using two band G-L theory was obtained in Ref. [16] with 
critical temperature parameters T c i=10 K and T c2 =3 K, approximately to yield a value of T c 
about 16 K. 

Note that the fitting between theory and experiement shows a deviation at 
temperatures very close to T c (see Fig. 2). It is likely that we may have two main reasons for 
the mismatch between theory and experiment at temperatures very close to T c . First, the 
thermal fluctuations may enhance the positive curvature near T c , which is ignored in the G-L 
theory as mean field theory. Second, in the solutions, we assumed a linear two band G-L 
theory for which there exists analytical solutions. For a better fitting, the cubic term must 
have been included in the two band G-L theory. Then the equations to be solve would 
necessiate the solutions for the vortex state in two band G-L theory. Our current interest is to 
seek solutions for vortex state and to map out the magnetic phase diagram in the vicinity of T c 
for the potential superconductor of MgB 2 . 

In summary, we have shown that the positive curvature of the upper critical field can 
be explained to a reasonable extent by seeking analytical solutions to the two-band G-L 
theory within an approximation. An exact fitting between experiement and the theory would 
necessitate exact solutions to the theory. However, the theory itself is phonemologic and it 
can explain the temperature dependence to some extent only at temperatures close to T c . 
Further work is therefore needed to elucidate the nature of superconductivity and the 
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interesting physical properties observed in MgB 2 at both microscopic level and the 
macroscopic level. 
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Figure Captions 



Figure 1. Temperature dependence of upper critical field H c2 (T). Open circle symbol shows 
the experimental data of Bud'ko et al. In Ref. [9], while filled circle symbols show the 
experimental data of Miiller et al. In Ref. [10]. The solid and the dashed-dotted lines are the 
best fitting curve and the two-band theoretical calculations of Ginzburg Landau theory, 
respectively. 

Figure 2. The first and second derivatives of Eqn.(10) and Eqn.(17) versus the reduced 
temperature. Open circle symbols show the best fit of the fitting formula of Eqn.(17), while 
the filled circle symbols show the derivatives by the G-L Eqn.(10). 
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Figure 1 
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